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We devise a theory of x-ray absorption by symmetric-top molecules which are aligned by an 
intense optical laser. Initially, the density matrix of the system is composed of the electronic ground 
state of the molecules and a thermal ensemble of rigid-rotor eigenstates. We formulate equations 
of motion of the two-color (laser plus x rays) rotational-electronic problem. The interaction with 
the laser is assumed to be nonresonant; it is described by an electric dipole polarizability tensor. 
X-ray absorption is approximated as a one-photon process. It is shown that the equations can be 
separated such that the interaction with the laser can be treated independently of the x rays. The 
laser-only density matrix is propagated numerically. After each time step, the x-ray absorption 
is calculated. We apply our theory to study adiabatic alignment of bromine molecules (Br2). The 
required dynamic polarizabilities are determined using the ab initio linear response methods coupled- 
cluster singles (CCS), second-order approximate coupled-cluster singles and doubles (CC2), and 
coupled-cluster singles and doubles (CCSD). For the description of x-ray absorption on the (jg Is ~+ 
(Tu 4p resonance, a parameter-free two- level model is used for the electronic structure of the molecules. 
Our theory opens up novel perspectives for the quantum control of x-ray radiation. 



PACS numbers: 33.20.Sn, 33.55.-b, 31.15.Qg, 33.20.Rm 
I. INTRODUCTION 

We study the interaction of an intense optical laser and 
hard x rays with molecules, a so-called two-color problem. 
The laser is assumed to be sufficiently intense to influ- 
ence molecular rotation, however, without vibrationally 
and electronically exciting the molecules in a gas sample. 
Specifically, we regard x-ray absorption by laser-aligned 
symmetric-top molecules. On the one hand, this leads 
to dynamical information of molecular rotation; on the 
other hand, this opens up novel ways to control x rays. 

Such a combination of an optical laser with x rays has 
proved very beneficial already in various setups. For in- 
tense laser pulses, ionization of atoms takes place, pro- 
ducing aligned, unoccupied atomic orbitals. They are 
probed by using XUV light or x rays to excite an inner- 
shell electron into them py, |2j, i^] . For slightly lower laser 
intensities, one only modifies the atomic Rydberg mani- 
fold leading to electromagnetically induced transparency 
for x rays 0, 0, @] . 

Laser-alignment of molecules was discovered exper- 
imentally by Normand et al. Ia| and was later an- 
alyzed theoretically by Friedrich and Herschbach Q- 
Other alignment techniques, however, did exist before, 
e.g., with strong static electric and magnetic fields 
for molecules with permanent electric and magnetic 
dipole moments, respectively [8|, |9|. Friedrich and Her- 
schbach [3] considered alignment with a continuous- wave 
laser which lead to eigenstates of the coupled problem of 
molecular rotation in the laser field. Ortigoso et al. [w\ 
carried out a systematic theoretical and computational 
study of laser-aligned linear molecules. They allowed for 
pulsed lasers and rotational wave packets to investigate 
the transition from a thermal ensemble to the continuous- 
wave solution of Friedrich and Herschbach Q- Aside of 
this so-called adiabatic alignment, there is also transient 



(or impulsive) alignment [8|, |9| which means that the 
molecules are "kicked" with an intense, short laser pulse 
and thereafter evolve freely. This leads to brief, recurring 
periods when all molecules in a sample are aligned [3, 0] . 
The work of Ortigoso et al. [l^] was extended recently to 
symmetric-top molecules by Hamilton et al. [ll|. 

The alignment of a gaseous sample causes the individ- 
ual molecules to basically rotate in the same way. Hence, 
at sufficiently low temperature, the x-ray absorption sig- 
nal reflects the properties of a single molecule rather 
than a rotational average. Adiabatic laser-alignment of 
molecules is the only choice when the molecules are to be 
studied with x rays from third-generation synchrotron 
sources due to long x-ray pulses of about 100 ps [I4I 
which are several orders of magnitude longer than the 
time scale of transient alignment of about Ips [1, @]. 
Particularly suited for the study of alignment is a de- 
pendence of the absorption cross section on the angle 
between the linear x-ray polarization and the flgure axis 
of symmetric-top molecules. Specifically, such a linear 
dichroism can be found on single, isolated resonances in 
the x-ray absorption near-edge fine structure (XANES) 
region [IJjImI- ^'^^ instance, above the carbon K edge, 
there is a a* shape resonance in hydrocarbons. The de- 
pendence of the cross section on the resonance on the 
orientation of the molecules was theoretically and exper- 
imentally studied with x rays for molecules adsorbed on 
surfaces |15{ . 

Using X rays to quantify molecular alignment has sev- 
eral advantages over conventional methods. Typically, 
molecular alignment is detected by, first, cleaving a 
molecular bond with one laser pulse and ionizing the 
fragments with a very intense other laser pulse [8]. The 
latter pulse leads to highly ionized fragments which un- 
dergo Coulomb explosion. Due to the complicated na- 
ture of the interaction of a laser with a molecule, this 



is a very intricate technique to study molecular align- 
ment; e.g., recently, the strong-field ionization proba- 
bility of aligned molecules was measured and a nontriv- 
ial, molecule-specific angular dependence was found [l6|. 
Therefore, this technique is prone to systematic errors 
and it is highly desirable to have an alternative method 
for the detection of molecular alignment. 

In order to theoretically investigate the two-color prob- 
lem, it is crucial to analyze and exploit the three different 
time scales of the involved physical processes which are 
assumed to be well-separated by a few orders of magni- 
tude from each other. First, the x-ray oscillations de- 
fine the fastest timescale. For example, at the K edge 
of a bromine atom, lox — 13474 eV [l3|, a cycle period 
amounts to — = 0.30694 as. Second, x-ray absorption 
leads to core-hole formation. A K vacancy in a bromine 
atom relaxes with a decay time of ^^ = 0.264 fs, i.e., the 
line width is F = 2.49 eV [17[. This relaxation defines 
the coherence time scale, i.e., the time scale on which x- 
ray amplitudes instead of intensities must be considered 
and thus interference effects in x-ray absorption may oc- 
cur due to molecular rotation. Third, the time scale of 
molecular rotation, and x-ray and laser envelopes; the 
quickest of them defines the third scale. This last scale 
varies with the choice of parameters; situations may arise, 
e.g., high laser intensity, where the separation of the third 
time scale from the second is no longer sufficient and a 
joint treatment of the processes on both time scales is re- 
quired. For this work, we expect the third time scale to 
be of a few tens of femtoseconds to several picoseconds, 
i.e., a separate treatment of the physical processes on the 
second and on third time scales is feasible. Because of 
this distinction of three time scales, the slower processes 
can be taken to be constant while approximations to the 
faster processes are obtained. Using the approximations 
to the faster process, the variation on the slower time 
scale can be determined. 

In this work, we study the influence of molecular align- 
ment on the absorption cross section for linearly polar- 
ized X rays. We devise a density matrix formalism for 
rotational motion of symmetric-top molecules under the 
influence of a laser. The laser-only equations are then 
extended to comprise the interaction with the x rays. 
Our work lays the theoretical foundation for our com- 
putational studies of bromotrifluoromethane (CFsBr) in 
Ref. y and the joint investigation of CFaBr with experi- 
mentalists in Rcf. 18. 

The paper is structured as follows. Section [Til discusses 
the theoretical foundation. We start with the laser-only 
problem in Sec. Ill Al the full two-color problem (laser 
plus x rays) is examined in Sec. IIIBI The x-ray absorp- 
tion cross section and derived quantities are discussed 
in Sec. Ill CI In Sec. IIIII a two-level electronic structure 
model is developed for bromine molecules. Computa- 
tional details are given in Sec. |IV] and the results are 
presented in Sec. El Conclusions are drawn in Sec. IVII 

Our equations are formulated in atomic units [l9|, 
where 1 hartree = 1 E^ is the unit of energy, 1 to is the 



unit of time, and 1 bohr = 1 oq is the unit of length. The 
Boltzmann constant fce is unity and lE^ ~ 3.15775 x 
10^ K is the unit of temperature. Intensities are given in 
units of lE/jfo^ag-^ = 6.43641 x lO^^Wcm^^. 



II. THEORY 

A. Laser alignment of symmetric-top molecules 

1. The freely rotating molecule 

We assume the Born-Oppenheimer approximation and 
the rigid-rotor approximation which implies that ro- 
tational and electronic degrees of freedom are decou- 
pled [1^, [201 and molecular vibrations are excluded [20] . 
Without interactions that reestablish a coupling, elec- 
tronic structure and molecular rotation can be treated 
independently. We deal in the following with two refer- 
ence frames: one is the laboratory or space-fixed frame 
where the coordinates are denoted by x, y, z. The other 
frame is the molecule-fixed or body-fixed frame with the 
coordinates a,b,c [20| . 

The molecular rotation is described by the rigid-rotor 
Hamiltonian which reads [201 



H, = A,T 



BJi 



CJ, 



c I 



(1) 



where A, B, and C are rotational constants. For a 
symmetric-top molecule, H^ commutes with the square 

of the angular momentum operator J^ , the z component 
of the angular momentum in the space- fixed frame Jz, 
and the c component of the angular momentum in the 
body- fixed frame Jc [20, [2l|, [22| . The three operators give 
rise to the rotational quantum numbers J, K , and M, re- 
spectively, with K, M = —J, . . . , J for integer J > 0. We 
assume throughout a prolate symmetric top, A > B ^ C, 
where the axis A is identified with the figure axis c. For 
an oblate symmetric top, A = B > C and the axis C 
becomes the figure axis c [20|. In this case, all A need 
to be replaced by C in what follows. We obtain from 
Eq. (|T]) the Hamiltonian 

H, = Bfi + {A- B)Jl (2) 

with the eigenvalues [20, [24] 

E.JK ^BJ{J+1) + {A- B) K^ . (3) 

The eigenfunctions of the time-independent Schrodinger 
equation with the Hamiltonian ([2]) are the symmetric- 
rotor functions [2^ 

(4) 
Here, (p, i?, x denote the three Euler angles that connect 
the space- and body-fixed reference frames [20, [22| and 



D'Lm -k^'-P^ ^1 x) represents elements of the Wigner ma- 
trix for angular momentum J [20, [2l|, [22| ■ The eigenfunc- 
tions of the general asymmetric rotor can be expanded 
in this basis @, Q . 

The molecules in a gas sample are initially not in 
a definite quantum state but are described by a ther- 
modynamic ensemble [20, l23|. Therefore, a density 
operator formalism is beneficial [23, [2J|. Let T be 
the rotational equilibrium temperature of the ensemble 
and /3 = ^. Then, the density operator of a freely rotat- 
ing symmetric-top molecule reads in its stationary eigen- 
basis dH) 



J,K,M 



J KM) wjK {JKM\ 



with the density matrix elements 



wjK ^ gi{J,K) 



-fSEjK 



(5) 



(6) 



The weights gj{J,K) reflect the nuclear spin statis- 
tics [23,[2| and 



J 



Z = ^(2J+1) J2 9iiJ,K)e-^^"< (7) 



J K=-J 

is the partition function. 



by q;^((jJl) in the space- fixed frame. 
Hamiltonian reads 28| 



1 



The interaction 



H^{t) = -^E^{t)^cx^{uj^)Edt) 



(10) 



After temporal averaging of Eq. (flOl) over an optical cy- 
cle, we obtain Hi^it) for the case cl = e^ as 



i/L(t) = -27ra/L(i)aLK) ■ 



(11) 



The dipole polarizability tensor is decomposed [2l|, |23 . 
|30| into a linear combination of a rank zero spherical 
tensor and the five components of a rank two spherical 
tensor [3l|; the zz component is rewritten as [30| 



Si 



K)=<oK) + <oK) 



(12) 



where we use a prime accent to indicate spherical tensor 
components. The term a^'Q{ujij) = iTra^(wL) is invari- 
ant under rotations and thus does not have any impact 
on molecular alignment. 

We would like to express a^^{uo-L) in terms of the 
dynamic polarizability tensor in the molecule-fixed 
frame ot^{uj-i) because the dynamic polarizability tensor 
is time independent in this frame. The transformation 
from the space-fixed to the body-fixed reference frame 
reads for spherical tensor operators [2l|, [23 | 



2. Interaction with a nonresonant laser 






(^l) 



E 



D^ 



,{^,d,x)a^:^,{LOi^) 



(13) 



Let the laser radiation of angular frequency wl be lin- 
early polarized along the z axis, i.e., the polarization vec- 
tor is gl = 6*2. The laser electric field is 



Eh{t) = ei,{t) CL cos(wL t) 



(8) 



Here, e-L{t) = ^/STToIiJt) is the envelope amplitude, a 
is the fine-structure constant, and Ihit) is the cycle- 
averaged intensity of the laser pulse at time t. We as- 
sume a Gaussian envelope with peak intensity /l,o and a 
[full width at half maximum (FWHM)] duration of tl 



T /j-\ T -41n2(-L; 



(9) 



The peak intensity of the pulse (at t — 0) depends on the 

X{0). 



In 2 -Epuia. 



laser pulse energy impulse as /l,o = '^ J ^ 
The factor X{0) = ^^^ represents the peak of a Gaus- 
sian radial profile of a FWHM width of pL- The treat- 
ment of X{g), with g being the radius perpendicular to 
the beam axis, would involve a number of computations 
with decreasing intensity. 

The interaction of molecules with an optical laser of 
a frequency far away from electronic transitions is de- 
scribed using the dynamic dipole polarizability of the 
molecules [20, [20, 123, 123 ! i* i^ given in terms of a sym- 
metric, rank two, Cartesian tensor [29l] which we denote 



We choose the axes of the molecule-fixed coordi- 
nate system to be the principal axes of a.^{uji,). 
Then, we have for a symmetric-top rotor, q;'^(wl) = 
diag[aj^(Li;L), Q;_L(a;L),a||(ti;L)]. There are only two dis- 
tinct diagonal elements: an (wl) is the polarizability par- 
allel to the figure axis of the molecule and a±{uji,) is 
the polarizability perpendicular to it. The components 
of the spherical tensors in Eq. p^ of rank zero and 
rank two 30] are found from the relations q;^q(ciJl) = 
i[2ai(c^L) + a||(c^L)], <o(^l) = |[a|| (c^l) - ^^(^l)], 
and a^±i(wL) = a^±2(^L) = 0. The following formu- 
las are expressed more compactly in terms of the average 
dynamical polarizability 



ai^h) 



It 



i—a.b,c 



i^(a;L) = ^a^(^L) + ^a||(^L) (14) 



and the dynamical polarizability anisotropy 

Aa{uji,) ^ a\\{uji,) - aj_{LUi,) . (15) 

Using Eqs. dnj, dUD, inD, HID, US]), and [13 



{JKM\DC„AJ'K'M')^^IJr^ 
xC{JJ"J'; KmK') C{JJ"J'; Mm' M') , 



(16) 



the matrix elements of the interaction Hamiltonian with 
res pec t to the symmetric-rotor functions ^ are found to 
be 



s pec t to t 
! [2ll,[2l 



{JKM\Hi^(t)\fK'M') 

4 
--na /L(i) Aq:(wl) Skk' hiM' 



2J+1 



C( J2 J'; KQK) C( J2 J'; MOAf ) 



(17) 



2J'-|-1_ 
27ra Ii^it) a(wL) 5jji Skk' Smm' 



Here, C{jiJ2J^]mim2m^) = (Jimi,J2m2 \j3m3} denote 
Clebsch-Gordan coefficients [2l|, [22|- The last summand 
in Eq. (fT7|l causes only a global energy shift; it has 
no impact on the rotational dynamics and may thus 
be dropped. Nota bene, the interaction matrix of a 
symmetric-top molecule for a linearly polarized laser is 
diagonal in K and M. 



Due to Eqs. ([3]) and (J18bp . the density operator in the 
Schrodinger picture simply results by multiplying the 
right-hand side of Eq. (|22p with the picture transforma- 
tion factor e~^'^-''J" * where [Eq. ^] 

coj'j" = Ej.K - Ej„K = B[J'{J' + 1) - J" {J" + 1)] , 

(23) 
for all K, i.e., the density matrix elements are 

9%^''\t) = {J'KM\g^Ut)\J"KM) 



Q -i'^j'j" * / f 



J'KM\Q^^,u{t)\J"KM) 



(24) 

The equation of motion of an ensemble of molecules in a 
laser field is given by the Liouville equation [23, |2J] 



i^nucl.l(i) = [-ffL,l(i),^nuclj(i)] • 



(25) 



3. Equation of motion for the rotational dynamics 

The time-evolution operator in the interaction 
picture — indicated by the subscript I — from tg to t is de- 
noted by ?7nuci,i(i, ts) [32] ■ The relation of a Schrodinger- 
picture state vector \il),t) and operator 0{t) to interac- 
tion picture state vector \il>,t)^ and operator Oi{t) is [2^ 



|^,i)i = c'^"* |V,t) , 
di{t) = c'^"' 0(t) e-'^"* 



(18a) 
(18b) 



with the Hamiltonian of the noninteracting system TJq- 
Here, -&o = H,.] later, from Sec. [HBJ onwards, H^ will also 
comprise the nuclear repulsion plus electronic Hamilto- 
nian JJj] of the molecule i^M, e.g., in Eq. (|35l) . 

At t = tsi there is no laser and the density operator 
of the ensemble of symmetric-top molecules in the inter- 
action picture results from Eqs. ((5|) and (J18b[) . For later 
times, the density operator is formally given by [2J| 

^nuclj(i) = t/nuclj(i, ^s) ^nuclj(is) ^nuclj(^' ^s) ■ (19) 

The time evolution can be recast as follows 



1 V-JifM, t)i = C/„ucu(i> is) I J KM) 



Let 



-J' J 



{t) = {j'KM\i>jKM.t\ 



(20) 



(21) 



be the time-dependent expansion coefficients. As the ma- 
trix representation (|17p of the interaction Hamiltonian 
with the laser is diagonal with respect to K and Af , so 
is the coefficient matrix ^. With Eqs. (O, (©, (fig)) . 
(|20)) . and (|2T|) . the density operator of a thermodynamic 
ensemble of symmetric-top molecules in a laser field has 
the following structure in the interaction picture: 



ftmcLI 



[t) - 



E 

J,J',J",K,M 
(KM) 



(KM), 



J'KM)c'j'^j">{t)wjK 



(22) 



Using Eqs. (|T8|) and (|22|) . one reduces Eq. ([25)1 to equa- 
tions of motion for the coefficients (1211) 



(26) 



ic^jff\t) - Y.{J'KM\Hi^{t)\J"KM) 



xe— *c(^f(0 



These equations form a system of coupled, first-order 
ordinary differential equations [33|. In order to solve 

Eq. (f26|) . the time integral f Cj,j {t') At' needs to be 

u 
found. By comparing Eq. (O with Eq. (|^^ . we obtain 

the initial conditions c}jij [ts) = 5jij with the help of 
Eqs. (O, (HOI), and (dH). Employing standard numerical 
methods for systems of differential equations [33| , an ap- 
proximate solution can be determined. In doing so, the 
full propagation interval is decomposed into A^t subin- 
tervals which are sufficiently short such that the coeffi- 
cients Cjnj (t) on the right-hand side of Eq. (|26p and the 
laser envelope, Eq. ^, can be assumed to be constant 
on each subinterval. 

To quantify the degree of alignment along the laser 
polarization axis, one typically uses the expectation 
value 8, 9|,[2l,[2|] 



(cos2i?)(t)=Tr[^nuci(i) cos^t?] 



(27) 



Here, ^ refers to the Euler angle that determines the tilt 
between the c axis and the z axis of the two coordinate 
systems [20, ^^- Using Eq. ([16]), the necessary matrix 



xc)j,/r*{t) {J"KM\ 



elements of cos^ z? = | Z)oo('^'''^' x) + ^ i^i the basis of 
symmetric-rotor functions ([3|) read 



( JKM\co^d I J'K'M' ) = 5kk' 5 mm' [\ ^j.r 
I J H^\ C{J2.r; KQK) C( J2 J'; Af OAf )] . 

(28) 



B. X-ray probe 

1. Interaction with x rays 

The x-ray pulse has the same functional form as the 
laser pulse ([8|) — however with L replaced by X — with the 
same, Gaussian, form for the envelope ex{t) as the laser 
pulse ([9]). The semiclassical Hamiltonian for the interac- 
tion of a molecule with x rays of angular frequency lux 
reads [si] 

Hx{t) = ~{d- ex.M)ex(t) cos(wxt) = hx{t) cos(cJx • 

(29) 

Here, d is the electric dipole operator. The x rays are 
assumed to be polarized linearly along the direction ex,M- 
Both vectors are given with respect to the molecule-fixed, 
Cartesian coordinate system. 

For simplicity, we transform the vectors from the 
Cartesian to the spherical basis [2l|, [22|- This change 
of coordinate system is indicated by a prime accent. It 
leaves the scalar product formally the same with the 
sum over Cartesian vector components j G {x, y, z} and 
a € {a, 6, c} replaced by a sum over projection quan- 
tum numbers m G { — 1, 0,1}. Using a coordinate system 
transformation [2l|, [22] analogous to Eq. p^ 



(^x.m) 



X.MM 



E 



Dt 



,(^,i9,x)(el). 



(30) 



we express the body-fixed x-ray polarization vector Cx m 
in terms of the space-fixed x-ray polarization vector e^. 
The dependence of Hx{t) in Eq. (P5)) on electronic and 
rotational degrees of freedom constitutes a coupling of 
both. Thus, the interaction with the x rays needs to be 
evaluated in terms of the direct product basis [35| 



\JKM,i) = \JKM)(g>\i) 



(31) 



consisting of the symmetric-rotor states | J KM ) and the 
electronic states \i). Here, « = denotes the electronic 
ground state. We obtain [ 36| 



{JKM,i\Hx{t)\J'K'M',i') 

-^li/ ■ SjKMj'K'M') ex(i) cos{ujxt) 



(32) 



The electric dipole matrix elements in the body-fixed ref- 
erence frame are designated by 

4, = (^|1'U'); (33) 

the nuclear rotation factor is defined by 

SjKMj'K'M' = ( JKM I el,M I J'K'M' ) . (34) 

The matrix elements (IMl) are easily evaluated using 
Eqs. O, dD, and dSOD- 



2. The two-color problem 

The joint Hamiltonian for nuclear dynamics, electronic 
structure, and the interaction with laser and x rays reads 

H = [H, + H^it)] ® iel + inucl ® Hm + Hx{t) . (35) 

Here, Id is the unit operator in the space of electronic 
states and Inuci is the unit operator in the manifold 
of symmetric- rotor states ^. Further, H,^ is given in 
Eq. ©, HL{t) in Eq. ^, and Hx{t) in Eq. ^. Finally, 
Hm denotes the fixed-nuclei molecular Hamiltonian [l^ . 
To describe the time-evolution of the interaction of a 
molecule with a laser and x rays, we extend the density 
operator formalism of Sec. lH A"3] to include electronic de- 
grees of freedom. Initially, the molecule is in the elec- 
tronic ground state |0). The corresponding stationary 
density operator simply reads [24| 



&i = |0)(0| 



(36) 



The initial rotational density operator is given by a 
thermal distribution ([5]). Then, the joint rotational- 
electronic problem is initially described by the direct 

product m, mi 



g == ^nucl ® ^cl 



(37) 



because, by assumption, for i < tg, we have Hx{t) — 
in Eq. ([35)1 . The relevant equations from Sec. Ill Al — 
Eqs. (USD, dSni), (ED), (USD, (BSD, and dMD— can now be 
generalized to comprise also the electronic degree of free- 
dom. To this end, we replace the nuclear rotational 
states I JKM ) with the direct product states | JKM, i ) 
from Eq. (j3ip [35|. Formally, this replacement means 
only that another index, namely, i needs to be added 
to the existing equations. The only thing one has to 
pay attention to is the fact that the interaction with 
the X rays — unlike the interaction with the laser ^Tf} — is 
generally not diagonal in K and M. Hence, the coeffi- 
cients (PT|) have to be replaced by 

cjKM,i;j'K'M',o(t) ^ {JKM,i\il;j'K'M',o,t)-^ . (38) 

By assumption, the molecule is initially in the elec- 
tronic ground state (i — 0). The energies of rotational- 
electronic states are 



E 



JK.i 



= E 



IK 



E.--lil 



<5o.) , 



(39) 



with EjK denoting the energy of nuclear rotational 
states ([U and Ei being the energy of electronic states 
of Hy[ in Eq. ((35|) . We consider x-ray absorption to pro- 
duce only electronically excited states (i ^ 0) with a 
single core hole which is an excellent approximation for 
the moderate x-ray intensities at third-generation syn- 
chrotron sources. The core vacancies relax by x-ray fluo- 
rescence and Auger decay [l^, [3 ■ The associated decay 
width is indicated as F; it is assumed not to vary for 



different core-excited states. This is an excellent approx- 
imation when the core vacancy is always in the same shell 
and has always the same angular momentum. 

We obtain the equation of motion for the joint problem 
in analogy to Eq. (^5)1 

i C,j'K'M'A';JKM,o{t) 

Y, {J'K'M\i'\Hi^{t)®i,, 

J",K",M",r" (40) 

+ Hx{t) I J"K"M", i")e'^J'K'y;j"K",^" * 

X Cj»K"M".i";JKMfi{'t) ■ 

In contrast to Eq. (P^ . here the sum extends also 
over K" , M" , and i" . The picture-transformation fac- 
tor [compare with Eq. ([23| ] becomes using Eq. ([39|) : 

^J'K',i':J"K",i" = Ej'K'A' — Ej'iK".i" 

r 

= Ej'K' — Ej^'K" + Ei' — Ei» — i— {5Qi>f — 5oi>) . 

(41) 



3. Recovering the laser-only equations 

Let us consider first the case i' = in Eq. (|^(I| . From 
Eqs. dni) and dSB) follows 



to the laser-only case ([26| , that the envelope of the x-ray 
pulse hx{t') on the right-hand side of Eq. P^ is approx- 
imately constant over the subintervals. (Please refer to 
the introduction concerning a discussion of the involved 
physical time scales which need to be taken into account.) 
For our forthcoming approximations to be appropriate, 
the duration of the subintervals are also required to be 
much longer than core-hole lifetimes i and the duration 
of x-ray oscillations — 

For molecules without permanent electric dipole mo- 
ment, i.e., dg'o — [Eq. p3|) ]. the second term in Eq. ((43l) 
vanishes exactly because the matrix element (j32p is zero. 
Otherwise, this term leads to a negligible integral on 
the time subintervals due to rapid oscillations, Eq. (j32p . 
of e"^-''^'.o;.7"R'".o * cos(a;x t). Thereby, we assume moder- 
ate x-ray intensities, i.e., the total x-ray absorption prob- 
ability for the full propagation interval is tiny. This also 
means that the third term is tiny because it depends 
on coefficients with i ^ 0. Omitting the third term in 
Eq. (P5|) implies that two- (and multi)-photon processes 
for X rays are neglected. 

With these arguments, the second and third terms in 
Eq. (P5|) are omitted leading to the laser-only equations 
of motion (j26p . In other words, we treat the interaction 
with the laser independently of the x rays in what follows. 



; J KM, i I Hi^it) I J'K'M', i' ) = 6k K' Smm' 4 
X {JKM\Hi^{t)\J'KM) 



(42) 



We arrive at 



i C,j'K'M',0;JKMfi{t) 

= Y, ( J'K'M' I HL{t) I f'K'M' ) 
J" 

X e'^-^'J" * C,jnK'M-fi;.JKMfi{t) 

+ Y ( J'K'M', I Hx{t) I J"K"M",0) 

J",K",M" 
X e'"'^'^'.o^J"^".o* Cj"A:-M",0;JKAf,o(i) 

+ J2 ( J'K'M', 1 Hx{t) I J"K"M", i" ) 

j",K",M",i" 
X e''^■'''<'.'>-■J"'<".^"UJ»K"M",^"■JKM.a{'t) ■ 

(43) 
The first term on the left-hand side of Eq. ([^5]) constitutes 
the interaction of the laser with the rotational degrees of 
freedom. It is the same as in Eq. (j26p in the expanded 
direct product space (PTjl . Similarly, the second term in 
Eq. (|43p represents transitions between rotational states 
due to the x rays. The third term in Eq. (|15)) describes an 
x-ray induced recombination of an excited electron with 
the core hole. 

The time evolution of this expanded system of 
differential equations (|15)) is treated in analogy to 

the laser-only equations in Sec. Ill A 31 by approximat- 

t 
ing J cj'K'M',0:JKM.o{t')dt' . For this purpose, we decom- 

pose [ts,i\ into A't subintervals. We require, in addition 



4- Electronic excitation due to x-ray absorption 

Let us now consider the case z' 7^ in Eq. (PD|) . With 
Eq. (I12D, we find 



( J'K'M', i' I Hx{t) I J"K"M",0) 



''iCj'K'M'.i';,JKM,o{t) 

= E 

J",K",M" 

+ J2 i J'K'M' I Hi^it) I J"K'M' ) 
J." 

X e'"^'''" * Cj"K'M'^t';JKM.o{t) 

+ Y ( J'K'M', i' \ Hx{t) I J"K"M", i" ) 

J",K",M",i" 
X e'"^''<''''--'"'<"''"*Cj'>K"M",t":JKM.o{t) ■ 

The meaning of the individual summands on the right- 
hand side of the equation is as follows. The first term rep- 
resents x-ray absorption by a ground-state electron lead- 
ing to electronic excitation of the molecule. The second 
term quantifies laser-induced transitions between rota- 
tional sublevels of an electronically excited state. Finally, 
the third term describes transitions by x-ray absorption 
between rotational-electronic core-hole states. Since we 
neglect x-ray nonlinear processes, we set this term to 
zero. One expects the second term in Eq. (|^^ to make 
a vanishingly small contribution, too, which is shown in 
the ensuing paragraphs. 

We need an approximate expression 

for cj" K' M' .i'\JKMfl{t) in the second summand on 



the right-hand side of Eq. (|44|) . We carry out the time 
integration of the right-hand side — thereby neglecting 
the third term — in analogy to the previous Sec. Ill B 31 
by decomposing the interval [ts,t\ into A't subintervals, 
letting t — tpf^ : 

C.J'K'M'.i'-.JKM,o{t) 

J",K",M" n=0 

X Cj/iK"M",Q:JKMfi{tn) ^ J' K' .i':,J" K" .o(^n,tn+l) 
-, Wt-1 

+ jJ2J2 ^ J'^'^^' I ^L(in) I J"K'M' ) 
J" n=a 

X Cj"K'M',i';JKM,o{tn) I j, j„{tn,tn+l) ■ 

(45) 

The remaining time integrals I j>j^, i'-j"K" oi^n, ^n+i) and 

Ij,j„(tn-,tn+i) in Eq. (|45|) can be carried out imme- 
diately; rapidly oscillating contributions are discarded 
employing the rotating- wave approximation (RWA) [3J| . 
We find 

^J>K',i'-J"K",oy^n,tn+l) 



cos 



{LUxt')e''^-'''<'-^'--'"'<"-o*' dt' 



2i {^,J'K'.i';J"K",0 — Wx) 

gU'^J'A",i';.7"R'",0^'^x) tn+1 

2i {i^.J'K',i'-.J"K",0 — ^x) 



and 



(46) 



Due to e-'^J'^'.''*4?^j„(t„,i„+i) w 0, the second 
term in Eq. (|45p vanishes in the Schrodinger picture 
because, by assumption, core-excited molecules decay 
rapidly on a subinterval [r(t„4.i — i„) 3> 1] and the laser- 
induced rotational response is negligible during the decay 
time of Y- 

Thus, Eq. (|45p goes over into 

Cj'if'Af ,i';J_ffM,o(i) 

« Y^ {J'K'M',i'\hx{t)\J"K"M",0) 



J",K",M" 



X Cj"K"M",0;JKM,oit) 



gi {l^jiK>,i';J"K",0-'^^) * 

-2 (wj'_ff',i';J"i<'",0 — Wx) 



(50) 
which is the final form for the coefficients describing elec- 
tronic excitations. We arrive at the corresponding ap- 
proximation of Eq. (j44p : 



Cj'K'M'S';JKM,o{t) 

i J2 {J'K'M\i'\Hx{t)\J"K"M\Q) (51) 

J",K",M" 
X e''^-'''^'-'^'"^".o*Cj"K"M",0;JKM,o(t) , 



which becomes identical to the time derivative of Eq. ([& 
after applying the RWA. 

In analogy to Eqs. (P^ and ([M|) . the Schrodinger- 
picture density matrix of the joint problem of molecular 
rotation and electronic excitation reads 



r(2) 



Ij'J"{tn,tn+l) 



e'^-"-'"'"{tn+l-tn) 



3iK/',7"t'dt' 



(47) 



In the last expression (I47p . the approxima- 
tion ujjijii (i,i+i — tn) ^ 1 is used, i.e., the rotational 
time scale is much longer than the duration of the small 
intervals. 

The coefhcients (pS]) are defined with respect to the in- 
teraction picture p^ . Yet the density operator eventu- 
ally will be used in the Schrodinger picture [see Eq. (|5^ ]. 
i.e., the coefficients involved are 

Cj'K'M',i';JKMfi{t) — e~' J'K'.i' Cj'K'M'.i';JKM.o{t) , 

(48) 
observing Eqs. ([M)) and (^1]) . The impact of this picture 
transformation factor on Eq. (I45p is exclusively given by 

its impact on the integrals Ij'K'i'-j"K"oi^ri^tn+i) and 

Ij,j„{tn,tn+i)- The first integral ((i5|) becomes with t — 

tNt- 

•= ^J'K'.i';J"K",Q\''ny''n+l) 

e-i(Sj/'K",o+'^x)t (49) 



(JjVt 



,n+l o: 



2i (wj'_ft:',i';j"i<-",o — wx) 



QjKM,i:J'K'M\i'{t) 

~ 2_^ Wj"K" CjKMS;J"K"M"fi{t) (52) 

J'\K"M" 
X Cj,xiM'i';J"K"M",o{t) ■ 



Explicit expressions for the density matrix elements are 
determined from Eq. (|52p by inserting the approxima- 
tion (j50p . One arrives at density matrix elements be- 
tween electronically excited states {i,i' ^ 0) 



QjKM,i;J'K'M'.i'{t) 

{J'K'M'.i' I /ix(i) I Jai^iMi, 0)* 



X 



LO 



.J'K',i'-,J2Kifi 



LUX 



(53) 

Observe that Qj j (i) from Eq. ([M]) enters here which 
contains the whole time evolution of the molecular ro- 
tation. Hence, the laser-alignment problem and the in- 
teraction with the X rays can be treated independently. 
The off-diagonal elements with either i = or i' = of 
Eq. ([5^ are all negligible due to rapid oscillations. 
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C. Photoabsorption cross section 

The probability of being in a core-excited state at 
time t follows from Eq. (|5^ and is given by 



Qit) — J2 Q.JKM,i;JKM,i{t) 
J.K.M 

— J2 f^J"K" \CjKA4,i;J"K"A4".o{t)\ : 

J.K.M, 
J" ,K".M" 

(54) 
The rate of change of the probability of being in a core- 
excited state then is 



with £y^{t) — y^Stt a lux Jx (t) and /x(i) = wx Jxit) [see 
Eq. dH)]. The same expression ([58|l is obtained from the 
theory of Ref. [I] by assuming that the laser-only density 
matrix is constant over short time intervals. Formally 
Eq. fSS]) resembles the expression (40) in Ref. ^ for the 
cross section of laser-dressed atoms. 

Given the instantaneous cross section ([58]) . we can 
compute experimentally more accessible quantities. The 
cross correlation of the cross section and the x-ray pulse 
is defined by [33| 



P^{r) 



<J^{t)Jx{t 



dt 



,th 



(59) 



g(t) = -rg(t) + rabs(0 



(55) 



in which 



abs 



(i) 



2 E 

J" ,K",M" 
X C* 



wj"K" Re [ e' 



-rt 



-JKM,i-J"K"M",oi''^) CjKM,t;J"K"M",a{t)] ■ 

(56) 
To interpret Eq. ([55]) . we assume, for the moment, no de- 
cay of core-excited states, i.e., F = 0. Then, the first term 
on the right-hand side vanishes and Q{t) = Fabs(i) holds. 
Hence, Q{t) quantifies the rate of change of the excited 
state population at time t due to the interaction with the 
X rays. In fact this means x-ray absorption because the 
nonlinear recombination process of the excited electron 
with the core vacancy is neglected in our equations [third 
term on the right-hand side of Eq. ([33])]. Consequently, 
TabslO is the rate of x-ray absorption at time t. With 
core- hole decay, i.e., F ^^ 0, Eq. ([55)1 has the form of a 



rate equation [3j|. On the right-hand side of Eq. ([55)) . 
the first term denotes the rate of excited-state popula- 
tion loss due to core-hole decay, whereas the second term 
represents the rate of increase of excited-state population 
due to x-ray absorption. 

From the x-ray absorption rate (|56p. the photoabsorp- 
tion cross section follows via 



a{t) 



abs 



it) 



Jx{t) 



(57) 



Here, Jx{t) = y^ represents the flux generated by A'x x- 
ray photons in the volume V that propagate in the 
same direction with the speed of light — [20| . The fine- 
structure constant is designated by a. We insert Eq. ([SO]) 
and the complex conjugate of Eq. ([?T|l into Eq. ([TO]) and 
use Eqs. (f24|) and ([57|) to obtain the instantaneous x-ray 
absorption cross section of a molecule in a laser field 



a{t) = 47rawxlin 



E 



AKM) 



it) E 



J,J',K,M 



J'KM,.J"K"M 



")idl^. 



J".K".M" 



^J"K"M",.JKM) 



E 



J"K'' 



Ee 



iT/2-EjK-Eo~ujx 



(58) 



Here, || and _L denote parallel and perpendicular laser 
and x-ray polarization vectors whereas th refers to a ther- 
mal ensemble, i.e., no laser radiation is present. Further, 
Piir) represents the total probability of x-ray absorp- 
tion for a time delay of r between laser and x-ray pulses. 
We define following ratios of cross correlations ([55]) for 
parallel and perpendicular polarizations and a thermal 
ensemble: 



R\\,^ir) 



^||/th(-r) 



R 



±/th 



(r) = 



P\\ir) 
P±ir) ' 
P\\{r) 
^th(r) ' 
P^r) 
Pthir) ■ 



(60a) 
(60b) 
(60c) 



The ratio (|60ap has been studied experimentally [18j . 

One can rewrite the cross section ([55)) to clearly exhibit 
its dependence on the angle between laser and x-ray po- 
larizations i^LX- We assume that laser polarization bl 
and x-ray polarization ex sue initially parallel and di- 
rected along the z axis. Keeping the laser polarization 
fixed bl =6*2, one obtains the rotated x-ray polarization 
vector v2(J\ with a passive rotation by ^lx around the 
y axis as 



ex = sin i9lx e^ -I- cos i^lx e^ 



(61) 



The vector ex needs to be transformed to the spherical 
basis [2l|, I24I and then be subjected to Eq. ([50)) to obtain 
the x-ray polarization vector in the molecular frame e^ m- 
Inserting this vector into Eq. ()34p . the angular depen- 
dence of the cross section ([58t reads 



(^i" ■ Sj'KM,J"K"M")id'j."0 ■ ^J"K"M".JKm) 
~ Sm WlX («0i" ' ^,J'KM,J"K"M")i"'i"0 ' ^J"K"M".JKm) 
+ cos UlX (rfoi" ' ^J'KM,J"K"M")i"'i"G' ^.1"K"M",JKm) 

+ sini?LX cosz?LX 

X [(dgj// • ■Sj'_R-Af^j"A'"A/")("i"0 ' ^J"K"M",JKm) 
+ (ripj/, • Sj'ifM,J"_R'"Af")("i"0 ' ^J"K"M",JKm)\ 

(62) 
with the assignments isjiKM,j"K"M")-m. = 

■^{J'KM\ - D\^ + Dl,^\J"K"M") and 



a) 



b) 



^ 



FIG. 1: (Color online) The two orbitals included in the two- 
level model for the electronic structure of a bromine molecule 
from which the x-ray absorption cross section of the molecule 
is determined [38| ]. a) Totally symmetric linear combina- 
tion (Tg Is (S^ symmetry) of the two Br Is orbitals; b) the 
(Tu 4p antibonding molecular orbital of E„ symmetry. 



(4'kM.7"K"M" )- = ( J' KM I Z?i „ I J"K"M" ) 

for m e { — 1,0,1}. The angular dependence of the 
cross section is expressed compactly by inserting 
Eq. §^ into Eq. dSH]): 



CT(t?LX,t) = sin2^LXCr_L(t)+COs2??LXCri|W 

-|-sin?9LX cos??Lxo-_L i|(i) , 



(63) 



where a±{t), a\\ (t), and ctj^ || (t) are the cross sections for 
parallel x-ray and laser polarization vectors, perpendic- 
ular polarizations, and a cross term. Generally, o'j_u(t) 
does not vanish. However, we will show that there are 
no cross terms for the special case of the two level model 
in the ensuing Sec. IIIII 

The angular dependence of a{'dj_,x,t) can be investi- 
gated experimentally [3|,|l8| in terms of a wave plate scan, 
i.e., the ratio 

oo 
-Rwavc(l?Lx) = -^^ / f7(l?LX,t)Jx(i)di 

— OO 

= cos^ t?LX i?||/th(0) + sin^ i?Lx i?±/th(0) 
+ sin?9LX cos?9LX^_L,||/th(0) 

(64) 
is measured, varying i?lx- In addition to Eq. (|60p . we 

define Rj_ n/th{r) = ^p-"/V using Eq. (ISi) with cr_L II (0- 



III. ELECTRONIC STRUCTURE MODEL 

We assume a two-level model for the electronic struc- 
ture of a bromine molecule (Br2). For simplicity, we treat 
the two states as Hartree-Fock orbitals [l9[ ; yet our rea- 
soning is general and not limited by this assumption. The 
model is justified by the experimental K edge absorption 
spectrum of Br2 which is shown in Fig. 1 of Ref. 4^; it 
exhibits an isolated pre-edge resonance which is assigned 
to a transition from the gerade linear combination of the 
Br Is atomic orbitals, the cTg Is molecular orbital, to the 
a^ Ap antibonding molecular orbital. The two molecular 
orbitals are displayed in Fig. [1] of this paper. 

We need to find out what the dipole selection rules are 
for the transition cTg Is -^ CTu 4p using symmetry argu- 
ments [4ll |43 . Bromine forms linear molecules of D^oh 



symmetry [42!]. The electronic transition is described 

by the matrix element dtj^is,a^ip = (cg Is | d |(Tu 4p) 
in Eq. (|33p formed with the Cartesian dipole opera- 
tor d = (a, 6, c)^ in the molecule-fixed frame. The van- 
ishing integrals rule [4l|, |42| says that the components 
of da is.(T„ 4p can only be nonzero, if the direct product 

representation of the representations to which (ffg Is |, d, 
and |(Tu4p) belong, contains the totally symmetric rep- 
resentation. The question therefore is, for i G {a^h,c\^ 
does r* '"s Is (g) F' (8) F'"" '^p contain F^^ ? 

There are two linear combinations of the Br Is orbitals. 
These are of EJ" and S^ symmetry, respectively; the 
CTu 4p orbital has S+ symmetry. The Cartesian compo- 
nent c of the dipole operator is of S^ symmetry whereas 
a and h are of Yiu symmetry. With the help of the direct 
product table of Dooh |i2l, we find that only the direct 
product of the combination of the Br Is orbitals with E+ 

symmetry with the c component of d and a^ Ap leads to 
a non- vanishing transition matrix element. This linear 
combination is thus identified with the (Xg Is molecular 
orbital in the two-level model. Only the c component dc 
oi da is,cru4p is nonvanishing. This implies that, accord- 
ing to the scalar product in Eq. ([32]), x rays are only 
absorbed with an x-ray polarization vector component 
along the internuclear axis; there is no absorption for 
perpendicular polarizations. In computations, we focus 
on ratios such as Eq. ([SO)) . Hence, our results are inde- 
pendent of the value of dc- Since the molecular orbitals 
are glued to the rigid nuclear framework, they rotate with 
the molecule. Hence, the photoabsorption is sensitive to 
the spatial alignment of the molecules. 

For the two-level model, one can further simplify the 
dependence on the angle between laser and x-ray polar- 
izations ^LX of the cross section ([S5]) . Due to the fact 
that the dipole vector d^ is,a-^ 4p has only a nonvanish- 
ing c component (or m = component) dc, the cross 
term in Eq. ([5^ becomes 



\dc\'^ sini^LX cos^ 

/ I 
+ s 



f.± 



LX [Sc,,J'KM,J"K"M" '=c,J"K"M",JKM 
c.J'KM,J"K"M" ^c.J"K"M",JKM\ ' 

(65) 



letting s 



',-L 
c,J'KM,J"K"M' 



\^J'KM.J"K"M")0 



and 



^c,J'KM,J"K"M 



^^ - (^^ 



J'KKLJ"K"M" 

j'KM,j"K"M")o- Recalling the 



rule P^ to evaluate the matrix elements in Eq. (|65p . 
we obtain conditions that need to be fulfilled such that 
the Clebsch-Gordan coefficients do not vanish. As we fo- 
cus on the m — component, we obtain from Eq. ([1 
the equality K — K" . Further, for s 
to be nonzero, M = 



'c,J'KM,J"K"M 



c.,J'KM,J"K"M" 

M" ± 1 must hold whereas 
„ ^ the assertion M ~ M" is nec- 



for 

essary. Both conditions cannot be fulfilled simultane- 
ously such that the cross term in Eq. (|65p and, there- 
fore, o-'^'"(i) in Eq. (p5|) vanishes. Our simplifications 
carry over to the signal of a wave plate scan; in Eq. 
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^L [Eh] 


Q [a.u.] 


Aa [a.u.] 


Q [a.u.] 


CC9. 


a [a.u.] 


crsD 


Basis set 


Aa [a.u.] 


Aa [a.u.] 


cc-pVDZ 





30.87 


37.76 


28.03 


33.91 


27.51 


32.40 




0.057 


31.28 


38.82 


28.43 


34.86 


27.89 


33.31 


aug-cc-pVDZ 





45.00 


30.99 


43.03 


27.85 


42.31 


26.47 




0.057 


45.68 


31.89 


43.79 


28.60 


43.04 


27.19 


aug-cc-pVTZ 





47.31 


30.66 


46.03 


27.57 


44.96 


26.41 




0.057 


48.06 


31.52 


46.87 


28.30 


45.76 


27.11 


aug-cc-pVQZ 





47.50 


30.72 


46.13 


27.74 


44.98 


26.73 




0.057 


48.26 


31.58 


46.99 


28.48 


45.79 


27.45 



TABLE I: Average dynamical polarizabilities q(cjl) and dynamical polarizability anisotropics Aq(cjl) of a bromine molecule 
in the light of a laser with photon energy uji^. The polarizabilities are determined with coupled-cluster linear response methods 
for several basis sets. Theoretical SDQ-MP4 values, Qrcf(O) = 44.79 a.u. and Aarcf(O) = 25.68 a.u., were taken from Table 3 in 
Ref . H 



-R_L,||/th(0) — and thus 

-Rwavo,modcl(l?Lx) = COS^ ^LX -R|| /th(0)+sin^ ^LX i?J_/th (0) ■ 

(66) 
with Eq. §U\). 

A two-level electronic structure model, similar to the 
one used in this section, is suitable when the x-ray en- 
ergy lies on an isolated resonance in the absorption spec- 
trum of a molecule. Neglecting the coupling between nu- 
clear and electronic motion as well as vibrations due to 
electronic excitation |44l ]. as we have done in this paper, 
the two-level model can easily be generalized to a full 
ab initio description because, in this case, the electronic 
processes are completely decoupled from the nuclear mo- 
tion. The electronic transitions matrix elements are de- 
termined in the fixed-nuclei approximation and inserted 
into Eq. 



IV. COMPUTATIONAL DETAILS 

Based on the theory of Sees. [TTl and IIIIl we have devel- 
oped the program ALIGNMOL of the fella package [43 ■ 
The program handles symmetric-top molecules where 
the special case of a linear molecule is treated by let- 
ting K = in all equations of this paper and eliminating 
sums over K. A large number of parameters are needed 
to specify the molecular properties, the numerical accu- 
racy, and the laser- and x-ray pulses which we discuss in 
the following paragraphs. 

The Br2 molecule is a linear rotor with an internu- 
clear distance of i? = 2.289 A [40|. The rotational con- 
stants in Eq. ^ are given hy A = and B = j^, 
where /i is the reduced mass [20]. There are two sta- 
ble, naturally occurring isotopes of bromine: ''^Br and 
^^Br [ii]. The first one, ™Br, has a mass of 78.9183361 u 
and a natural abundance of p 793^. = 50.69% whereas 
the second one, *^Br, has a mass of 80.916289 u and a 
natural abundance of psiBr = 49.31% [H, \^. This 



leads to three possible combinations of isotopes in Br2 
and thus to three different rotational constants: B — 
0.081537 cm-i for ™Br2, B = 0.079524 cm^^ for ^iRra, 
and B = 0.080530 cm-i for ^Br-^^Br. 

Both bromine isotopes have a nuclear spin quantum 
number of / = | [40] ■ The statistical weights in 
Eq. ^ for an isotope-pure diatomic molecule with nu- 
clear spin /, are given by [20, l2a] 



gi{J,K) = 



1 ; J even 
; J odd 



(67) 



For ^^Br-^^Br, the nuclei are distinguishable and the sta- 
tistical weight is gi{J,K) = 1. As the equations in this 
paper are mostly nonlinear in the rotational constants 
and the statistical weights, it is required that for each 
of the three combinations of bromine isotopes, a sepa- 
rate computation of the desired quantities is carried out. 
The results are then averaged assigning the weights p ra^^ 
to results for ™Br2, p\i^^ to ®^Br2, and 2p79Q,.p8iQ,. 
to 79Br~8iBr. 

We assume a rotational temperature of Br2 of T = 
30 K; the initial rotational density matrix (0 needs to 
be represented in terms of a finite number of symmetric- 
rotor functions (|4]) [for K — because Br2 is linear]. 
We chose angular momenta up to Jmax = 70. The 
value for Jmax is determined to be sufficiently large such 
that the contribution with J — Jmax to the partition 
function ([7]) is < 10~^. For each J — 0, . . . , Jmax all 
symmetric-rotor functions ^ with M — ~J,...,J are 
employed to represent the density operator; sums over K 
are omitted. The total energy of the Gaussian laser 
pulse (HI) is £;puiso = 2mJ; it lasts (FWHM) tl = 100 ps. 
The energy and duration of the pulse imply a peak in- 
tensity of fa 10^^ -^ on the axis of a laser beam with 
a Gaussian radial profile perpendicular to the beam axis 
of a FWHM width of ql = 40 /xm. We do not consider a 
radial intensity variation but take the value on axis. 

As the molecule is in the field of the aligning laser, 
we need to determine the dynamical polarizability tensor 
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FIG. 2: (Color online) The time evolution of (cos^ ■&) (t) 
[Eq. p7|) ] of bromine molecules is given by the solid (red) 
line as a function of the laser pulse which is depicted by a 
dashed (black) line. 
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FIG. 3: (Color online) Ratio of cross correlations R^y±{T) 
[Eq. (|60a|l ] for bromine molecules as a function of the time 
delay between laser and x-ray pulses. 



in the molecule-fixed reference frame a (wl) [27|, |48|. 
This means specifically the average dynamical polar- 
izability a{LUh) [Eq. (|14[) ] and the dynamical polariz- 
ability anisotropy Aa{uji^) [Eq. (J15p ] which both enter 
Eq. (|17|1 . The quantities can be computed in fixed- nuclei 
approximation utilizing the ab initio quantum chemistry 
program package DALTON [49|. We employ molecular 
linear-response theory in conjunction with wave func- 
tion models of increasing sophistication: coupled-cluster 
singles (CCS), second-order approximate coupled-cluster 
singles and doubles (CC2), and coupled-cluster singles 
and doubles (CCSD) [4g|. The cc-pVDZ [correlation con- 
sistent polarized valence double zeta] and aug-cc-pVXZ 
[augmented cc-pVXZ] (for X = D^, Q [D: double, 
T: triple, Q: quadruple]) basis sets [50, [5l| are used to 
describe the bromine atoms. We examine the static case, 
wl — OE/i, and the dynamic case for the photon energy 
of an 800 nm Ti:Sapphire laser, ojl = 0.057 E/j. 

Our average polarizabilities and the polarizability 
anisotropics are listed in Table ID A moderately rapid 
convergence with respect to the basis set quality is 
observed for the individual methods; the difference is 
much smaller between corresponding results for aug- 
cc-pVTZ and aug-cc-pVQZ than between the results 
for ujl — OEft and o^l = 0.057 E/^, i.e., the basis set 
limit is nearly reached. We conclude that the remain- 
ing deviations between methods for a specific wl and 
the aug-cc-pVQZ basis reflect the differences between 
the coupled-cluster methods. Comparing our best results 
for wl = OEft,, CCSD and aug-cc-pVQZ, with theoreti- 
cal SDQ-MP4 values from Table 3 in Ref. [H, we note a 
satisfactory agreement. In our computations, we use the 
value of Aq;(0.057E;i) — 27.45 a. u. for the polarizability 
anisotropy. 

The coefficients (PT|) of the rotational density matrix 
in Eq. ([^ are represented here by iVcq = Jmax values 
for J' which determines the number of equations of mo- 
tion (p6)) for each initial symmetric-rotor state | J KM). 
We use the Runge-Kutta algorithm of fourth order [33| 
to compute the time evolution of the density matrix via 
Eq. (P5|) . The propagation takes place over the time in- 
terval [— 400ps,400ps] which is subdivided into A^timc — 
6002 time steps. 

A Gaussian shape of the x-ray pulses similar to the 



laser pulse ([9]) is assumed. As the interaction with 
the X rays is treated in terms of a one-photon process, 
the specification of a value for the peak intensity be- 
comes unnecessary; the relevant prefactor cancels once 
ratios are formed as in Eqs. 160]) and (|66l). Only, the 
FWHM duration of the x-ray pulse needs to be given. We 
choose Tx = 120 ps, a typical value for a third-generation 
synchrotron radiation source [12| . The cross correlation 
of the cross section and the x-ray flux ([59)1 is computed 
via numerical integration using the trapezoid rule [33| . 
To this end, we assume the laser pulse to be centered at 
zero, Eq. ([5]), and use an equidistant grid of iVx = 1003 
x-ray pulses, i.e., 1003 different time delays r, over the 
whole propagation interval [—400 ps, 400 ps] . 



V. RESULTS AND DISCUSSION 

Here, we exemplify the theory of the previous sections 
for bromine molecules. We picked laser and x-ray pulse 
parameters in the previous Sec. IIVI the pulse durations 
were chosen such that they are comparable; the resulting 
molecular dynamics is shown in Fig. [21 We give the time 
evolution of (cos^ d) (t) for the rotational wave packet 
[Eq. P7)) ]. From the figure it becomes clear that we are in 
the regime of adiabatic alignment because the temporal 
evolution of ^cos^-i^Wi) and of the laser pulse are prac- 
tically the same [^, [3 ■ Even the shape of the molecular 
response is very similar to the one of the laser pulse. Yet 
the response of the molecule is not totally Gaussian as 
indicated by the slight deviation between the two curves 
in Fig. [1 

The ratio of cross correlations -R||/j_(t) [Eq. (|60ap ] be- 
tween cross section and x-ray flux is displayed in Fig. [31 

The ratio of total x-ray absorption between laser on 
and laser off is given by Eq. ([66)1 . It is plotted in Fig. [4l 
in relation to the angle between laser and x-ray polariza- 
tion vectors i?lx- Without laser, there is only a thermal 
ensemble and the angular dependence vanishes. 



VI. CONCLUSION 

We devised a theory for the x-ray absorption of 
symmetric-top molecules subject to an intense nonres- 
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FIG. 4: (Color online) Dependence of the ratio between the 
x-ray absorption signals of laser- aligned bromine molecules 
versus a thermal ensemble on the angle i?lx between laser 
and x-ray polarization vectors. 



onant laser. The rigid-rotor approximation was assumed 
and a density-matrix formalism was used to describe a 
thermal ensemble of rotational states evolving in time. 
The x-ray absorption is treated in terms of a one-photon 
process. The very different time scales of molecular rota- 
tion and electronic relaxation (core-hole lifetime) facili- 
tated a separation of the two problems, i.e., the molecular 
rotation was treated separately from the x-ray absorp- 
tion. The electronic structure of bromine molecules was 
described in terms of a simple parameter-free two-level 
model for the tJg Is -^ Cu 4p pre-edge resonance. The 
symmetry of the resonance caused a linear dichroism — 
light was only absorbed for a linear polarization compo- 
nent along the internuclear axis of Br2 — which lead to a 
modulation of the photoabsorption with respect to the 
alignment of the molecule. 

Our theory was realized in terms of the ALIGNMOL pro- 
gram which is part of the fella package [4^. The dy- 
namic dipole polarizability for the laser frequency was 
determined for Br2 using three ab initio coupled-cluster 
methods of increasing precision: CCS, CC2, and CCSD; 
it was found not to deviate significantly from the static 
polarizability. We investigated the time evolution of the 
alignment of Br2 molecules in the adiabatic regime. The 
x-ray absorption is studied in terms of the experimen- 
tally accessible ratio of cross correlations between cross 
section and x-ray pulse for several combinations of laser 
on and off and parallel and perpendicular linear polar- 



ization vectors of the laser and the x rays. A continuous 
variation of the angle between the polarization vectors 
was also investigated. 

This article forms the theoretical basis for future 
studies. In forthcoming papers, we will apply our 
theory to the symmetric-top molecule bromotrifluo- 
romethane (CFsBr) [3, [lH. There, we will compare our 
results with experimental data. Further perspectives in- 
clude obtaining absolute cross sections from our theory. 
For this purpose one needs a quantitative description 
of the electronic structure of the molecule, i.e., one re- 
quires transition dipole moments instead of our two-level 
electronic structure model. The moments can be ob- 
tained from existing ab initio program packages. Our 
work provides an alternative way to the established but 
intricate [iq Coulomb-explosion technique ^ to study 
laser-controlled molecular rotation in the gas phase which 
additionally opens up the possibility to study the laser- 
manipulated rotation of an ensemble of molecules in so- 
lution |52l | where the Coulomb-explosion technique is 
bound to fail. Moreover, with our theory, one can in- 
terpret experimental data; the symmetry of molecular 
resonances in the gas phase enters our formulas in terms 
of selection rules for dipole matrix elements. In conjunc- 
tion with experiments, our approach offers a novel route 
to determine the symmetry of resonances. It thus com- 
plements conventional methods using clamped molecules 
on surfaces [53| and angle-resolved photoion yield spec- 
troscopy |54| . 
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